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Abstract 

^~~i The theory predicts that the spin-wave hfetime tl and the hnewidth of ferromagnetic resonance 

^N /S.B can be governed by random fields and spatial memory. To that aim the effective field around 

, ^ which the magnetic moments perform a precession is superimposed by a stochastic time dependent 

'^ magnetic field with finite correlation time. The magnetization dynamics is altered by inclusion 

1-^ of a spatial memory effect monitoring a non-local interaction of size ^. The underlying Landau- 

,^ Lifshitz-Gilbert equation (LLG) is modified accordingly. The stochastic LLG is equivalent to a 

Q Fokker-Planck equation which enables to calculate the mean values of the magnetization vector. 



S 



Within the spin-wave approximation we present an analytical solution for the excitation energy 
and its damping. The lifetime and the hnewidth are analyzed depending on the strength of the 
— ^ random field D and its correlation time Tc as well as the retardation strength Fq and the size ^. 

Q Whereas tl decreases with increasing D, retardation strength Fq and t^, the lifetime is enhanced 

O 

'— ' for growing width ^ of the spatial retardation kernel. In the same manner we calculate the exper- 

►^ imentally measurable hnewidth A.B is increased strongly when the correlation time Tc ranges in 

^r^ the nanosecond interval. 
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I. INTRODUCTION 

Ferromagnetic resonance (FMR) is a powerful technique to study magnetic materials, in 
particular the inherent magnetization dynamics [U [2] . So the observable FMR-linewidth is 
very sensitive to the underlying dynamical processes as well as the real structure of the ma- 
terial like anisotropy. From a theoretical point of view the Landau-Lifshitz-Gilbert equation 
(LLG) [21 H], see Eq. ([I]) in the present paper, is an appropriate tool to investigate magnetic 
excitations and dissipative processes as the damping of the excitations. Although, the LLG is 
known since a few decades it is still a standard model to analyze magnetodynamics. Recently 
the Gilbert damping parameter was determined experimentally for ferromagnetic thin films 
in [S] and by first-principle calculations for itinerant ferromagnets in |B]. Obviously, the ap- 
plicability of the LLG depends on the physical situation in mind. In case the magnetization 
is not conserved the Landau-Lifshitz-Bloch (LLB) equations are more appropriate, in partic- 
ular in the vicinity of the phase transition as demonstrated in |7|. The LLB equations were 
used to investigate magnetization switching near the Curie temperature in [8l|9]. Moreover, 
the geometrical configuration of the sample are able to play an important role in measuring 
the FMR-linewidth. Related to this fact the contribution of the Gilbert damping to the 
linewidth can be superimposed by extrinsic effects as magnon-magnon scattering processes 
[To] which become of the same order of magnitude or even exceed the Gilbert damping. Es- 
pecially for an in-plane configuration where the external field as well as the magnetization lie 
in the film plane the influence of two magnon processes to the FMR-linewidth cannot be ne- 
glected im [12] . Those theoretical results predicting a nonlinear dependence of the linewidth 
on the frequency were extended to the case when the magnetization is tipped out of plane 
[13j . Different experimental findings emphasize the importance of extrinsic contributions for 
in-plane setups, see [HHTT]. A quantitative separation of Gilbert damping and two magnon 
scattering contributions was carried out [lll[T6l[T7j. Contrary to these observations there are 
other investigations y^, which offer no qualitative difference between in-plane and normal- 
to-plane measurements. In both realizations the linewidth depends linearly on the frequency 
even for frequencies smaller than 10 GHz . Such theoretical and experimental works suggest, 
among others, that the FMR-linewidth is exclusively controlled by the Gilbert damping and 
exhibits a pure linear frequency dependence in a perpendicular configuration with respect 
to thin films measurements. Furthermore, the two magnon scattering is supposed to be of 



less importance in bulk ferromagnets [16]. Thus the LLG equation seems still applicable 
to describe magnetization dynamics provided the physical situation is carefully analyzed as 
pointed out in [2] . A more realistic magnetization dynamics requires a modification of the 
LLG. Recently, the anisotropic damping and its manifestation in the FMR-linewidth has 
been discussed by several authors [T9H2T] . An alternative formulation of Gilbert damping 
by means of scattering theory was discussed in [22]. In addition, ferromagnetic resonance 
measurements were used as well to investigate spin transport in magnetic single and double 
layer structures [23] • Moreover, very recently it was shown that the transfer of spin angular 
momentum can induce ferromagnetic resonance dynamics in a ferromagnetic film due to 
the spin Hall effect in an adjacent film with strong spin-orbit scattering [21]. Related to 
this phenomena it was reported on the direct time-resolved measurement of spin torque in 
magnetic tunnel junctions to detect resonant magnetic precession due to an oscillating spin 
torque [25] . A theory of ferromagnetic resonance in perpendicular magnetized nanodisks is 
suggested in pS]. 

To push forward the theory stochastic forces and non-local interactions should be included 
into the model to gain a more realistic description of magnetic materials and to reveal unex- 
pected behavior as for example the noise suppression by noise behavior argued in [27]. The 
effects of noise in magnetic nanosystems obeying spin torque dynamics are investigated in 
[28] [29] . Experimentally, the role of noise in magnetic systems was prospected in [301 E] ■ 
The present work is addressed to the influence of randomness on the magnetization dynam- 
ics. As the two new aspects the system considered is simultaneously subjected to feedback 
coupling and to a stochastic fleld with colored noise. The starting point is the LLG equation 
which is generalized in a manner that both spatial memory effects and a temporal stochastic 
fleld with a flnite correlation time is incorporated into the model. Previously the influence 
of colored noise p^ and retardation effects [33] within the LLG were analyzed separately. 
Otherwise, both effects can occur simultaneously. Consequently we study a combined model 
concerning both kind of impacts, feedback and randomness. As demonstrated in former 
papers there exits the possibility that the total damping, originated by the Gilbert damp- 
ing and that one induced by memory effects are able to cancel by the distinct damping 
mechanisms. In this paper we are interested in the FMR-linewidth. The corresponding 
parameters range in such reasonable intervals where different dissipation sources are not 
observable. The main goal is to calculate the FMR-linewidth and to discuss its dependence 



on the parameters characterizing randomness and retardation. 

Let us give a brief outhne of the paper. In Sec. 2 we present the mathematical model and 
its underlying basic assumptions. The stochastic LLG equation is equivalent to a Fokker- 
Planck equation which is derived approximately in Sec. 3. This equation enables to compute 
the mean values of the magnetization. The results are discussed in detail in Sec. 4. Finally, 
we conclude by summarizing the results and by an outlook in Sec. 5. 

II. MODEL 

As already indicated in the introduction we are interested in micro- and nanosized mag- 
nets. Therefore a coarse-grained description is an appropriate tool to investigate magnetic 
material. In this mesoscopic description the discrete magnetic moments are replaced by 
a spatiotemporal vector field m(r,t). The interaction and the dynamics of the moments 
are formulated in a continuous approximation. The situation is schematically illustrated 
in Fig. [I| Here, the magnetization m(ri) represents the magnetic properties within the 
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FIG. 1. Illustration of the of the coarse-grained mesoscopic model. The Si represent microscopic 
magnetic moments which are related to the magnetization m. Further explanation can be found 
in the text. 



volume-element d^ri which is build around the position ri. The field m(ri) stands for the 
total set of microscopic spins which will be visible if one zooms into the microscopic struc- 



ture. The huge number of microscopic degrees of freedom within d^ri are substituted by 
a single degree of freedom, namely the mesoscopic quantity m which can be considered as 
the sum over the microscopic spins located at equivalent crystal positions. Moreover, the 
magnetization vector field m(ri) is assumed to be oriented continuously in space. The basics 
of our model consists in this mesoscopic description discussed before. Further, the system is 
supposed to offer an uniaxial anisotropy where the direction of the anisotropy axis is denoted 
by I/. Our calculations refer to weak excitations which evolve as spin waves and possess a 
finite life time. Both quantities are found in the long wave-length limit ga ^ 1, where q is 
the amount of the wave vector and a is the lattice constant. This assumption reflects the 
mesoscopic level of description. Experimentally the dynamic behavior of the magnetization 
m can be detected for instance by means of ferromagnetic resonance (FMR). Because the 
main goal of the paper is to put forward the modeling towards more realistic systems we 
develop a dynamic model for the magnetization field m(r) in which retardation effects as 
well as stochastic fields are included. In particular, the aim is to relate our findings for the 
magnetic excitations and their damping to an experimentally accessible quantity, namely 



the FMR-line width AB, cf. Eq. (23). 



As underlying model we start from is the Landau-Lifshitz-Gilbert equation 



dm 7 

dt 1 + a2 



Beflf + a [m X Beff] , (1) 



which will be generalized accordingly. In Eq. ([I]) the quantities 7 and a are the gyromagnetic 
ratio and the dimensionless Gilbert damping parameter, respectively. In this description 
m(r, t) is the unit vector m = M/M^ with the magnetization M the saturation magneti- 
zation. The local effective field Beff(r, t) causes the precession of the magnetization. In 
general, the effective field Beg is composed of different contributions, an internal field due 
to the interaction of the spins, the magnetic anisotropy and an external field. This effective 
field can be derived from the Hamiltonian of the system by functional variation with respect 
to m 



B.,^-M-|^. (2) 



The Hainiltonian T-i can be expressed as 



"H = j (fr {Wex + Wan + Wext] 



Wex = -Mj{Vin)\ (3) 

Wan = -MsK Sin^ 9 , Wext = "Bext " M . 

The quantities J = J a? and K designate the exchange energy density and the magneto- 
crystaUine anisotropy energy density. Here, J is the couphng strength between nearest 
neighbors referring to the isotropic Heisenberg model [3S] and a is the lattice constant. 
Further, Bgxt is the static external magnetic field. The quantity 9 is the angle between 
the the local magnetization m and the anisotropy axis u = (0,0,1). We assume that u 
points in the direction of the easy axis in the ground state. Therefore, i^ > characterizes 
the strength of the anisotropy. In deriving Eq. ^ we have used m^ = 1. Let us stress 
again that this assumption seems to be correct if the temperature is well below the Curie 
temperature [7]. Our calculations based on the LLG suggest that other damping mechanism 
such as an extrinsic magnon-magnon scattering due to magnetic inhomogeneities should be 
inactive and hence they are irrelevant. In thin films this situation can be achieved when 
both the magnetization and the static external field are perpendicular to the film plane. 
In our model this situation is realized when both the easy axis of the anisotropy u as well 
as the external field Bgxt point into the z-direction. Hence the equilibrium magnetization 
is likewise oriented parallel to the z-axis. This situation corresponds to a normal-to-plane 
configuration. From here we conclude that the application of the LLG leads to reasonable 
results. For a different realization an alternative dynamical approach seems to be more 
accurate, see also the conclusions. 

To proceed further, the vector m is decomposed into a static and a dynamic part termed 
as i-i and i/' = (V"!; "^2; V's); respectively. In the frame of spin wave approximation we make 
the ansatz 

m(r, t) = n + il}{r,t) = fiv + il^ , fi = const. , (4) 

Combining Eqs. (|2| and ^ yields the effective field 

Beff = JV'V'-^V'' + Bext , V' = (V'1,^2,0). (5) 



It is appropriate to introduce dimensionless quantities: 

IB I ^^^ 

K 

The quantity Iq is called the characteristic magnetic length [36] whereas the parameter 
e reflects the ration between the strengths of the external and the anisotropy field. For 
convenience later we will substitute t — > t again. So far we have introduced the LLG in 
Eq.Q in its conventional form and incorporated our special basic model assumptions for 
a ferromagnetic material below its Curie temperature. To proceed toward a more realistic 
description of magnets the LLG will be extended by the inclusion of retardation effects 
and random magnetic fields. Whereas retardation is implemented by a memory kernel 
r(r,r';t,t') a stochastic field r]{r,t) contributes additionally to the effective field, i.e. 

Beff(r, t) -^ beff(r, t) = Beff(r, t) + r]{r, t) . (7) 

Taking both effects into account we propose the following generalized LLG 

^^^ = fdt'Jd'r'Tir-r';t-t') 

x(--^m(r',t')x [Ksir',t')+ (8) 

[ 1 + a^ 

+a[m(r',t')xbeff(r',t')]] 
where the stochastic field is included in the dimensionless effective field as 

beff = /oV'V-V'' + ^bo + r7(r,t). (9) 

The unit vector bo indicates the direction of the external magnetic field. In general, the 
kernel should respect the retardation concerning temporal and spatial processes. More 
precise, a change of the magnetic moment at position r should influence another moment at 
position r' and vice versa. This influence is thought to be an additional contribution which 
should not be confused with parts of the exchange interaction in the effective field, i.g. 
the length ^ on which spatial retardation effects are relevant could be of a different order of 
magnitude in comparison with the lattice constant a. Insofar, a purely coordinate dependent 
part of the kernel reflects a kind of non-local interaction. All moments within a radius i^ 



contribute to the interaction. Likewise a temporal feedback mechanism can be taken into 
account due to the fact that the transport of information from one magnetic moment to its 
neighbors needs at least a finite albeit small time. Such an in-time retardation mechanism 
is considered already in [33j- Here we concentrate on instantaneous retardation in time 
whereas the spatial part is realized for simplicity by a Gaussian shape 



r(r;t) = 5(t)|^-^exp 



, 2 

r 



(10) 



where Fq and ^ determine the strength and the size of the retardation, respectively. The 6- 
function in the last equation signalizes that all contribution to the interaction within a sphere 
with radius ^ contribute simultaneously to the interaction. As discussed below a typical value 
for ^ is assumed to be of the order 10~®m, i.e. the time for the signal propagation within ^ 
is about 10^^^ — 10~^^s. Because this time is much smaller as the lifetime of the spin-waves, 
see the discussion below, we conclude that delay effects within the region with radius ^ can 
be neglected. As indicated in Eq. ([t]) the noise ri{r,t) can also depend on space and time, 
i.e. in general random forces can effect the value of the magnetization at different positions 
in a distinct manner while additionally their fluctuations are also time dependent. Such 
a behavior maybe lead back to local infinitesimal temperature gradients or defects. The 
random field ^(r,t) is regarded as a colored noise the statistical properties of which obey 
the following relations 



Dal3 



T, 



exp 



0/3 



t-t' 



T, 



a/3 



(11) 



T"q/3->0 



>2D^p5{t-t') 



The components ria{t) have a zero mean and a finite correlation time. As an aside in 
the limit r — )■ the usual white noise properties are recovered. However, we want to 



concentrate on the more realistic colored noise case with r > 0. In Eq. (11) we assume 
ria{r,t) = ria{t). In other words the total system is affected by the same random influences. 
This may be reasonable if we have a well controllable constant temperature over the whole 
sample and an ideal sample without defects. Let us briefly summarize the new properties 
of the model defined by Eqs. ([s]) and ^. The influences of retardation and a multiplicative 

8 



noise as well are implemented in the conventional Landau-Lifsliitz-Gilbert equation. After 
a general incorporation into the model we had to limit the properties of both retardation 
and noise to an idealized situation in order to obtain analytical results in the subsequent 



section. However, although each of the Eqs. (10) and (11) represents a simplified version 



of a more general case the linking between both by means of the equation of motion for 
the magnetization in Eq. ([s]) models a quite complex behavior which is partly indicated in 
Fig. [2] While the exchange interaction is a short range coupling over a lattice constant a, 



oc J 



feedback oc Tr 




t+n 



FIG. 2. Schematic depiction of the difference between the exchange interaction J and the coupUng 
due to retardation oc Fq. As is visible feedback mechanisms can range over a larger distance 
S, ~ P = na, where n is integer. 



the interplay due to retardation with strength Fq can cover a distance p which is a multiple 
of the lattice constant. If this distance p is comparable to the characteristic length scale ^ in 



Eq. (10) retardation effects should be relevant. This microscopic picture can be transferred 



to a mesoscopic one and means a kind of non-local interaction. On the one hand at every 
spatial point the same kind of noise affects the magnetization. Otherwise, the magnetization 
m(r, t) takes different values at distinct positions r and therefore, the impact of the noise 
might be slightly different, too. Although spatial alterations of the noise are not regarded in 



the correlation function defined in Eq. (11) the memory kernel respects spatial correlations 



within ^ as seen in Eq. (10). Insofar the effect of noise at different spatial positions is 



transmitted by the memory kernel F(r, r'; t,t'). Another important hallmark is that the 



noise-noise correlation function Xap{t^t') is featured by a finite lifetime, cf. Eq. (11). For 
the forthcoming calculations we assume that r„^ = Tc^ap- Likewise the matrix of the noise 



correlation strength is supposed to be diagonal, i.e. Dki = DSm- Hence, the two important 
stochastic parameters are the correlation time r^ and the correlation strength D whereas 
the relevant parameters originated by the retardation are the retardation strength Fq and 



the retardation length ^, see Eq. (10). The results will be discussed in terms of the set of 
parameters D,Tc,To and ^. 



III. STATISTICAL TREATMENT 

Eqs. ([s]) and ^ represents the stochastic LLC Due to the coupling to the stochastic field 
T7(r, t) the magnetization field ni(r, t) becomes a stochastic variable. To calculate the mean 
values of m one needs the probability distribution P(m, t). To that aim the current section 
is devoted to the derivation of an approximated Fokker-Planck equation which allows to 
find the equations of motion for averaged quantities. To that purpose let us reformulate the 
model presented in Eqs. Q, dsl) and ([9]). After Fourier transformation ^^(q, t) = J-'T{'?/'(r, t)} 
we find in linear spin-wave approximation 

— (^„(q, t) = Aa[(p{q, t)] + B^^[(p{q, t)] r7^(t) 

The vector A and the matrix B posses the components 

/ 

. fiq,0 



(12) 



a^ 



\ 



-(/3/i + e) {aj^ifi + 1P2] 

{f3fi + e) (</?! - afnp2) 





(13) 



/ 



and 



/. 



B 



l + a2 



A 



(14) 



-<-P3 a/iv?3 (pi-afnp2 

Here the function /(q, ^ is the Fourier transform of the memory kernel r(r,t) defined in 
Eq. (10) and /i and (3 are introduced in Eqs. (Ill) and ([6]), respectively. Notice that / depends 



/ 



only on the absolute value q of the wave vector and takes the form 

1 



/(?, = To exp 



2^2 



e^g 



(15) 



To get the probability distribution function of the stochastic process determined by Eqs. (11) 



and (12)-(14) we define according to [371 EH] 

Picp,t) = {6[cpit)-cp]) 



(16) 
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Here the symbol < ... > means the average over all realizations of the stochastic process. As 
usual ip{t) represents the stochastic process whereas (p are the possible realizations of the 
process at time t. Due to the colored noise the corresponding Fokker-Planck equation can be 
obtained only approximatively in lowest order of the correlation time. The time evolution 



of Eq. ( 16 ) can be written in the form 

d 



dt 



Picp,t)=CP{cp,t) 



(17) 



In deriving this expression we have used the time evolution of (p{t) according to Eq. (12), 
the Novikov theorem [39] and the correlation function given by Eq. (11) with Taj3 = Tc6af3, 
Dap = D6ai3- The form of the operator C is given in a correlation time and cumulant 
expansion while transient terms have been neglected [^0] - H2] 

r\ f^ f) 



X 



D[B^f,{ip)-TcM^p{ip 



D^Tr 



1 

2' 



d 



d 



(18) 



+ ^B^^M—K^,f}^{ip) 



with 






A, 



OB. 



7/3 



B 



dip 
dB. 



/t/3- 



7M 



dA, 

d(Pu 
dB^p 

-D, 



(19) 



fiiy • 



Notice that summation over double-indices is understood. The single probability distribution 



is determined by the operator C in Eq. (18) which enables us to find the equation of motion 



for the expectation values (v^a). It follows 

d 



dt 



{iPa{t))={Aa) + D 




(20) 



2 \dipy \ dip^ 
Notice that in the white noise case all terms oc r^ would vanish. 
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IV. RESULTS AND DISCUSSION 



We find an analytical solution for the colored noise problem in Eq. (20) by standard 



Greens' function technique and Laplace transformation. After performing the summation 



in Eq. (20) while making use of Eqs. (13), (14) and the expressions in Eq. (19) the result 
reads 



/ 



{^(t)) 



St , 



e~^^cos{nt) e-^*sm{nt) 



\ 



V 



• iV'o) 



(21) 



7 



-e-^^sm{nt) e"^*cos(fit) 
e-'^* 

where (<^o) = {^{t = 0)) are the initial conditions. Physically, the parameters 5, a and Vt 
play the roles of the inverse magnon lifetimes and the frequency of the spin wave, respectively. 
They are determined by 



5 =« /i [e + /3 /i] ^f^^ + D [2 - a V'] 



l + a' 



+ 2 D Tc a /i [£ + /3 /i] 



+ D Ve [1 - 6 a^ /i^] 



/Ml 

l + a2 



/Ml 

l + a2 

/(g,0 

l + a2 



fi 



[e + p n\ z;—^ — ^ + 3Da n ' 



1 + a' 



+ Dtc [a^ jj^ -l][e + /3 jj,] 



1 + a^ 



+ 



a =2D 



D^Tr^ 



2 ,,2l 



a/i[ll — 3a /i 



f{q,0 
l + a2 

f{q,0 
l + a2 



(22) 



fiq,0'- 



a^ 



[4 /:> Tc a /i (e + /3 /i)] 



f{q,0 '' 
l + a2 



+ D^ Te [3 a'jj' + 1] 



f{q,0' 



1 + ^2 ■ 

Note that the parameters of the retardation mechanism, the strength Fq and the length scale 



^, are included in the function f{q,^) defined in Eq. (15). The two important parameters 
originated from the noise are the correlation time r^ and the correlation strength D of the 



random force. Both affect the quantities in Eq. (22) as well. We proceed by studying the 



system under the variation of these four model parameters. To be comparable to FMR 
experiments we refer to the following quantities 



TL = {6-fK)-\ AB = 1.16 



au 



7 



l.lGaKQ 



(23) 
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i.e. the lifetime tl of the spin waves and the FMR-hnewidth AB, compare [H [T7], which 



are related to the dimensionless inverse lifetime 6 and frequency il from Eq. (22). Here the 
frequency u is tantamount to the resonance frequency of the spin waves. The lifetime tl 
and the linewidth AB are given in Sl-units. Notice that the frequency independent part 
ABq, typically added on the right-hand side of the equation for AB is already subtracted in 



Eq. (23). The contribution ABq is supposed to take into account magnetic inhomogeneities. 
For a quantitative evaluation we need to set the model parameters to reasonable values. In 
doing so we also refer to Eq. (|6]). First let us start with fixed values. For the Gilbert damping 
parameter we choose the bulk value for Co which was found to be a ~ 0.005 |131 HI]. A 
similar value {a ^ 0.0044) was measured for a FE4/V4 multilayer sample in perpendicular 
configuration where only intrinsic Gilbert damping is operative |16j . 

The anisotropy field K is estimated as follows. Since the exchange interaction is typically 
about 10^ times larger than relativistic interactions which are responsible for anisotropy 
[15] and the magnetic exchange field can adopt large values we estimate the anisotropy 
as i^ = 0.1 T. Since we are interested in small excitations transverse to the anisotropy 
axis ly we suppose fi = 0.9 for the time independent part of the magnetization pointing in 
the direction of the anisotropy axis u, compare Eq. (|4]). Moreover, the gyromagnetic ratio 
7 ^ 1.76 X 10^^ (Ts)~^. The characteristic magnetic length defined in Eq. ^ is of the order 
Iq ~ 10~^ m. For the calculations a static magnetic field of about 0.5 T is taken into account 
which corresponds to the scaled external field e = 5. Notice that the dispersion relation Q 



in Eq. (22) is g-dependent. In the following we assume a medial value q = 10^ m~^. The 
parameters which are altered in the upcoming analysis are the noise correlation strength D 
and the retardation strength Fq. We investigate our model for both values ranging in the 
interval [0,10]. After this estimation the two parameters, the noise correlation time Tc and 
the retardation length ^ are left over. For a comprehensive estimation we suggest that ^ is 
ranged in 10~^^m < ^ < 10~^m. The lower limit is smaller than a typical lattice constant 
a c^ 10^^*^ m where the upper limit is a few orders of magnitude larger than the lattice 
constant. Likewise the correlation time Tc captures a quite large interval. Remark that 
we keep the notation Tc, especially with regard to Fig. |3| although we also designated the 
dimensionless correlation time as Tc. In order to cover a wide range the time interval is chosen 
in between atto- and nanoseconds. The results are depicted in Fig. [3] and Fig. |4j In Fig. [3]the 
behavior of the FMR-linewidth AB as well as the lifetime of the spin waves tl, introduced 
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FIG. 3. The FMR-linewidth and the Hfetime depending on: (a) the noise correlation strength D for 
Tc = 568 as, To = 1, C = lO^'^m; (b) the noise correlation time Tc for D = 1, Tq = 1, ^ = lO^^m; 
(c) the retardation strength Fq for Tc = 568 as, D = 1, S, = 10^^ m; (d) the retardation length S, for 
Tc = 568 as, Fq = 1, Fq. The other parameters take Iq = 10^^ m, q = 10^ m~^, e = 5, fi = 0.9 and 
a = 0.005. 



in Eq. (23), are shown in dependence on the different model parameters explained above. 



The inffuence of the correlation noise strength on AB and tl is shown in Fig. [sta). Whereas 
the linewidth decreases only very weak linearly when the noise strength D is increased, 
the lifetime of the spin waves tl reveals a strong dependency on D. This is indicated by 
the fact that tl is monotonic decaying while it covers several orders of magnitude with 
growing noise strength D. The curve shape for the lifetime tl seems to be comprehensible 
because the stronger the stochastic forces are correlated and interact mutually the faster the 
coherent motion of the spin moments is destroyed. This microscopic picture is reasonable 
under the premise that the evolution of spin waves is based on the phase coherence between 
adjacent magnetic moments. Apparently the frequency and consequently the linewidth AB 



show only a quite small effect, compare Eq. (23). Therefore, the variation of D reveals 



no significant influence on the frequency velocity of the moments. A distinct behavior is 
depicted in Fig. [sFb) for AB and tl as a function of the noise correlation time Tc. Both 
the linewidth and the lifetime remain constant for large interval of the correlation time Tc, 



14 



roughly speaking for Tc ranging from as to ps. If the correlation time is in between ps and ns 
the linewidth AB increases about a factor of 20 and the lifetime r^ decreases to a value about 
9-times smaller. Thus t^ affects both AB and t^ in an opposite manner provided the noise- 
noise correlations occur on time scales larger than ps. In this regime a growing correlation 
time Tc implicates likewise an enhancement of the resonance frequency of the spin waves 



CO oc AB, see Eq. (23). Simultaneously the spin wave lifetime tl declines strongly. Such a 
behavior may be attributed to a 'stochastic acceleration' which on the one hand enhances 
the frequency but on the other hand drives neighboring magnetic moments out of phase 
coherence. Remark that for times Tc > Ins the linewidth AB tends to infinity. This effect is 
not shown in the picture. Concerning the influence of the retardation parameters we refer to 
Fig. IsFc), which illustrates the influence of the retardation strength Fq. As recognizable the 
FMR-linewidth exhibits a seemingly linear dependence as function of Fq while AB grows 
with increasing retardation strength. The lifetime t^ decreases in a non-linear manner. The 
decay covers a range of ~ 3 orders of magnitude. We suggest the following mechanism 
behind this effect: Let us consider two moments both localized at arbitrary positions within 
the retardation length ^ as schematically displayed in Fig. |2j The mutual coupling due to 
retardation between both characterized by Fq leads to a phase shift between neighboring 
spins. Therefore, the phase coherence originated by the self-organized internal magnetic field 
is interfered in view of an interplay within the feedback coupling in coordinate space. The 
stronger this interaction Fq is the faster is the damping of the spin waves. Accordingly, spin 
wave solutions for different values of Fq are plotted exemplary in Fig. |4j The retardation 
length ^ influences AB and ti as well as is visible in Fig. [sld). Here the quantities AB and 
Tl remain constant for a retardation strength ^ ranging within the pm regime and a few 
tenth fim. For larger ,^-values the linewidth AB decreases while the lifetime tl increases. In 
the regime ^ > l/xm the linewidth AB tends to zero and the lifetime r^ — ;■ oo. This behavior 
is not depicted in Fig. [std). Notice that for reasonable values of ^ which not exceed the 
sample size the mentioned situation is not realized. The shapes of the curves in Fig. p[d) 
may be explained as follows. This graph corresponds to a fixed retardation strength Fq 
while the retardation length ^ is enlarged. Again we refer to the physical picture where 
the internal field, originated by the mutual interaction of the moments, and the coupling 
due to the retardation operate as opposite mechanisms. The interplay happens in such a 
manner that an increasing retardation strength Fq weakens or destroys the phase coherence 
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FIG. 4. Evolution of spin waves for different values of the retardation strength Tq. The other 
parameters take Iq = 10~^ ui, q = 10^ m^^, e = 5, /i = 0.9 and a = 0.005, Tc = 568 as, D = 0.5 and 
^ = 10"^ m. 

between adjacent spins. Yet it is found that a growing retardation length ^ counteracts the 
damping of the spin waves. As a consequence we suppose that the more spins are involved 
into the retardation effect, i.e. the larger the parameter C, becomes, the more the damping 
is reduced. In other words it seems that retardation effects can average out if sufficiently 
many magnetic moments are involved. 



V. CONCLUSIONS 

In the present paper we have studied a model on a mesoscopic scale realized by means 
of Landau-Lifshitz-Gilbert dynamics. The magnetization is driven by an effective magnetic 
ffeld. This field consists of an internal field due to the exchange interaction, an anisotropy 
field and a static external field. Additionally, the effective field is supplemented by a time 
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depending random one obeying colored noise statistics. Moreover, the stochastic LLG is 
generahzed by the introduction of a retardation kernel depending on the spatial coordinates 
only. Such a kernel simulates a kind of non-local interaction of size ^. After deriving an 
approximated Fokker-Planck equation we were able to calculate the mean values of the 
components of the magnetization in the linear spin wave approach. They depend strongly 
on the parameters characterizing the retardation (strength Fq, length ^ ) as well as the 
stochastic (strength D, correlation time t^ processes. As a result of the analysis we found 
that the increase of the retardation strength Fq compared with the growth of the retardation 
length ^ can entail conflictive effects on the lifetime r^. The main results are depicted in 
Fig. |3] There, in addition to the lifetime of the spin waves tl the FMR-linewidth Ai? is 
displayed. In doing so we want to provide comparability to experimental investigations 
based on ferromagnetic resonance for cases when the LLG is applicable. Let us remark 
that also other mechanisms are able to contribute to the damping process. As suggested in 
[21 US] the Bloch-Bloembergen equations |171 HH] are more appropriate for in-plane configu- 
rations in thin films. These equations are characterized by two relaxation times. Another 
approach with different relaxation processes is based upon the Landau-Lifshitz-Bloch equa- 
tions IIHI [50] . Our method including the inevitable stochastic forces can be likewise applied 
to the modified set of equations. 

One of us (T.B.) is grateful to the Research Network 'Nanostructured Materials' , which 
is supported by the Saxony-Anhalt State, Germany. Furthermore, we are indebted to Dr. 
Khalil Zakieri (MPI of Microstructure Physics) for valuable discussions. 
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